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ABSTRACT 

We present a statistical study of substructure within a sample of ACDM clusters 
and galaxies simulated with up to 25 million particles. With thousands of subhalos 
per object we can accurately measure their spatial clustering and velocity distribution 
functions and compare these with observational data. The substructure properties of 
galactic halos closely resembles those of galaxy clusters with a small scatter in the mass 
and circular velocity functions. The velocity distribution function is non-Maxwellian 
and flat topped with a negative kurtosis of about -0.7. Within the virial radius the 
velocity bias b = cTsub/cDM ^ 1.12 ±0.04, increasing to b > 1.3 within the halo centers. 
Slow subhalos are much less common, due to physical disruption by gravitational tides 
early in the merging history. This leads to a spatially anti-biased subhalo distribution 
that is well fitted by a cored isothermal. Observations of cluster galaxies do not show 
such biases which we interpret as a limitation of pure dark matter simulations - we 
estimate that we are missing half of the halo population which has been destroyed 
by physical overmerging. High resolution hydrodynamical simulations are required to 
study these issues further. If CDM is correct then the cluster galaxies must survive the 
tidal field, perhaps due to baryonic inflow during elliptical galaxy formation. Spirals 
can never exist near the cluster centers and the elliptical galaxies there will have little 
remaining dark matter. This implies that the morphology-density relation is set before 
the cluster forms, rather than a subsequent transformation of disks to SO's by virtue 
of the cluster environment. 

Key words: methods: N-body simulations - methods: numerical - dark matter — 
galaxies: haloes — galaxies: clusters: general 



1 INTRODUCTION 

Early simulation work that attempted to follow the merg- 
ing hierarchy produced a final dark matter structure 
that was nearly entirely smooth IWhite 1976.: White et alJ 
| 1987 : Carlbere 1994; Sum mers. Davis fc EvrardI Il995t 
iTormen. Bouchet fc Whitelll997ir The reason for this be- 
haviour was debated in the lit erature as being due to physi- 
cal or numerical ov ermerging ( White & Roes 1978; "Carlberi 
Il994l: Ivan Kampen 1995; M oore. Katz fc Lake 1996). The 
development of fast algorithms to accurately integrate the 
orbits of millions of particles overcame this problem. The 
first halos simulated with sufficient resolution contained of 
the order a thousand substructure ha los with propertie s 
that resembled galaxies within clusters iMoore et al.lll99a) . 
These simulations took many months using parallel gravity 
codes running on hundreds of processors. 

Ongoing research in this area has given many interest- 
ing results and we list some of the main conclusions here 
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jGhiffna et alj Il998l 


lOkamoto & Habd |l998|: iKlvoin et al 


'l999a 


VKlvoin et al. 


1999b. Moore et al." 1999;'Ghiena et al 


.2000; 


Springel et al. 


2001; Do Lucia et aL 2004): (i) Subha- 



los make up a fraction of between 5 and 10% of the mass 
of virialised halos. (ii) Halos on all mass scales have similar 
substructure populations, (iii) The mass and circular veloc- 
ity function of subhalos are power laws with slopes -1 and 
-3. (iv) Velocity bias between the subhalos and smooth dark 
matter background may be significant, (v) The radial num- 
ber density profile of subhalos is shallower than the dark 
matter background, (vi) Subhalos are significantly rounder 
than field halos (vii) The orbits of subhalos are close to 
isotropic with apo:peri approximately 4:1. (viii) Subhalos 
suffer mass loss from tidal stripping which modifies their 
outer density profiles, (ix) The tidal radii of subhalos de- 
creases with cluster-centric position, (x) Most of the surviv- 
ing population of subhalos entered the parent halo late. 

Several of these statements remain controversial and 
further work is necessary to clarify certain issues. In this 
paper we re-address conclusions (i)-(v) and attempt to an- 
swer some of the remaining questions, including: What is 
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the scatter in the mass and circular velocity distributions? 
Is there a positive or negative v e locity bias and if so what is 
its ori gin? lOhiena et al.l l)200CJl : IColin. Klvpin. fc Kravtso" 



20001) claim a positive velocitv bias whilst ISpringe^^J 
20011) report a negative velocity bias. Have we converged 



in the properties of subhalos, including their radial dis- 
tribution and mass functions? The inner regions (r < 0.2 
Tviriai) of clusters aud galaxy dark matter simulations are 
nearly smooth but is numerical overmerging still occurring 
in these very high density regions? Does the spatial dis- 
tribution of galaxies in clusters resemble that of the sub- 
halos in simulations? On galaxy scales the observed dis- 
tribution of satellites is more concentrated than the sim- 
ulations. Theory can be reconciled with the observations 
if it is assumed that the vis ible satellites are a biased 
subset of the total population ([Taylor. Silk, fc BabuJl2003l : 
iKravtsov. Gnedin &: Klvpinll200jr On cluster scales we do 
not expect to find "dark galaxy halos" therefore it is inter- 
esting to compare the observed distribution of galaxies with 
the distribution of substructure. 

In this paper we analyse a sample of six high resolution 
simulations of clusters containing between 5 and 25 million 
particles integrated with high force accuracy. We compare 
the mass functions with a sample of galactic mass halos 
with slightly lower resolution. These new simulations are 
presented in Section 2 and the general properties of the sub- 
halos are given in Section 3. 



2 NUMERICAL EXPERIMENTS 

Table gives an overview of the simulations we present in 
this paper. With up to 25 x 10^ particles inside the virial 
radius of one cluster they are among the highest resolution 
ACDM simulations performed so far. They represent a ma- 
jor investment of computing time, the largest run was com- 
pleted in about 5 10^ CPU hours on the zBox supercomputer 



2.1 N-body code and numerical parameters 

The simulations were carried out using PKDGRAV w ritten 
by Joachim Stadel and Thomas Quinn l)StadeJl200ll ). Indi- 
vidual time steps are chosen for each particle proportional to 
the square root of the softening length over the acceleration, 
Ati — rj^yrjoi. We use rj — 0.2 for most runs, only in run 
D9lt we used larger timesteps rj = 0.3. The node-opening 
angle is set to 9 — 0.55 initially and after z = 2 to 6 = 0.7 
to speedup the runs. The code uses a spline softening length 
e, forces are completely Newtonian at 2e. In table^eo is the 
softening length at z = 0, emax is the maximal softening in 
comoving coordinates. In most runs the softening is constant 
in physical coordinates from 2; = 9 to present and constant 
in comoving coordinates before, i.e. emax = lOeo- In runs C9 
and Fdcm the softening is constant in comoving coordinates 
for the entire run, in run F9ft the softening has a constant 
physical length for the entire run. 
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2.2 Initial conditions and cosmological parameters 

We use a ACDM cosmological model with parameters from 
the first y ear WMAP results: = 0.732, f7„ = 0.268, 
as = 0.9, (ISpergel et al J 120031) . The i nitial conditions ar e 
generated with the GRAFIC2 package iBertschingeill20o'lh . 
Six clusters were selected from a parent simulation and res- 
imulated with much higher mass and force resolution, de- 
tails about the selection and the refinement are given in 
iDiemand etai\ (l2004b|) . We label the six cluster (ordered 
by their mass) with letters A to F and with a number that 
gives the refinement factor in length relative to the 300'^ 
in (300 Mpc)'' parent simulation, e.g. 'D12' is the fourth 
most massive cluster in our sample, and the mass resolu- 
tion corresponds to (12 x 300)^ particles in a 300 Mpc cube 
simulation. 

We also present results from four medium resolution 
galaxy mass halos which we label GO, Gl, G2 and G3. These 
halos contain 2-4 million particles within their virial radii. 
The parent simulation is a 90 Mpc cube resolved with 300'^ 
particles initially. The four galaxies all lie within a volume 
of about 1000 cubic Mpc (at z = 0) which was refined by a 
factor of 12 in length to reach the resolution given in Table 

m 



2.3 Substructure Identification 

Within the virial radius of the high resolution CDM simu- 
lations we can resolve thousands of substructure halos, i.e. 
self-bound over-dense clusters of particles (See Figure 1). 
They span a wide range in mass, from the resolution limit of 
a few tens of particles up to few percent of the cluster mass, 
i.e. from 1O*M0 to 1O^''M0. Some of the subhalos even con- 
tain their own substructure. Therefore robust identification 
of subhalos a very difficult task, there is no general, param- 
eter free method that is able to extract the entire hierarchy 
of halos. 

\\ . identify subhalos with SKID dStadell I200I and 
with a new parall el adaptive Friends of Friends ('FoF', see 
iDavis et al.] ^8^) group finder ('AdFoF'). SKID calculates 
local densities using an SPH kernel, then particles are moved 
along the density gradient until they oscillate around a point 
(i.e. move less than some length Then they are linked to- 
gether using FoF with this I as a linking length. AdFoF first 
calculates the background density of the cluster pbg using 
spherical bins. The linking lengths for the particles are set 
to fe = (ApBc/n^p)"^''^, rnp is the particle mass, A = 5 is 
the density contrast, the only free parameter of this method. 
Two nearby particles can now have different linking lengths, 
they are considered as friends if one of them considers the 
other one his friend, i.e. the maximum of the two linking 
lengths is used. Both the SKID and the AdFoF groups are 
checked for self-boundness and unbound particles are re- 
moved with the same iterative procedure. 

We compared SKID results (using I = 1.5eo, I = 4eo 
and I = lOeo) with the AdFof results and we also visually 
compared the identified groups with the density map of the 
cluster: SKID with I = 4eo adequately identifies the small- 
est subhalos and the centers of the largest subhalos. For 
the latter the calculated bound mass is underestimated. Us- 
ing I = lOeo can be cure this, but then some of the small 
subhalos are missed. The AdFoF has the advantage that 
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Table 1. Parameters of resimulated clusters. The last four columns give properties of all subhalos with at 
least 32 bound particles, their number, bound mass fraction, the radius of the innermost subhalo and the 
velocity bias b = (Tsub/croM- In clusters A9 and C9 these structures are the cores of massive clusters that 
are about to merge with the main cluster at z = 0. 



Run 


[kpc] 




1O15[M0] 


virial 

[kpc] 


''^v,max 

[km s^] 


^vc.max 

[kpc] 


"halo 


Mvirial 


^sub,min 

[kpc] 


b 


^9 


2.4 


24'987'606 


1.3 X lO^^^ 


2850 


1428 


1853 


5114 


0.07 


126* 


1.10 


B9 


4.8 


11'400'727 


5.9 X 10^^ 


2166 


1120 


1321 


1940 


0.12 


162 


1.12 


C9 


2.4 


9'729'082 


5.0 X lO^"* 


2055 


1090 


904 


1576 


0.11 


77* 


1.15 


D3h 


1.8 


205'061 


2.8 X lO^** 


1704 


944 


834 


36 


0.03 


260 


1.11 


D6h 


1.8 


1'756'313 


3.1 X lO^** 


1743 


975 


784 


307 


0.04 


136 


1.11 


D6 


3.6 


1'776'849 


3.1 X lO^"* 


1749 


981 


840 


322 


0.05 


227 


1.13 


D9 


2.4 


6'046'638 


3.1 X lOi"* 


1752 


983 


876 


929 


0.06 


126 


1.11 


D9lt 


2.4 


6'036'701 


3.1 X 10^* 


1752 


984 


841 


912 


0.05 


183 


1.11 


D12 


1.8 


14'066'458 


3.1 X lO-"^"* 


1743 


958 


645 


1847 


0.06 


136 


1.11 


E9 


2.4 


5'005'907 


2.6 X 10^"' 


1647 


891 


889 


829 


0.06 


172 


1.11 


F9 


2.4 


4'567'075 


2.4 X 10" 


1598 


897 


655 


721 


0.06 


176 


1.08 


F9cm 


2.4 


4'566'800 


2.4 X 1014 


1598 


898 


655 


661 


0.06 


127 


1.08 


F9ft 


2.4 


4'593'407 


2.4 X lOi"* 


1601 


905 


464 


706 


0.06 


161 


1.07 


GO 


0.27 


1'725'907 


1.01 X 10^2 


260 


160 


52.2 


144 


0.03 


16 


1.05 


Gl 


0.27 


1'905'113 


1.12 X 10^2 


268 


162 


51.3 


189 


0.04 


20 


1.03 


G2 


0.27 


3'768'008 


2.21 X 10^2 


337 


190 


94.5 


462 


0.04 


21 


1.10 


G3 


0.27 


2'626'202 


1.54 X 10^2 


299 


180 


45.1 


314 


0.03 


28 


1.12 



in principle it links together all particles in regions with a 
density contrast of A against the background density. With 
A = 5 AdFof finds the same groups as SKID, but the current 
version using the spherically averaged density for the back- 
ground also finds some spurious groups since the background 
isodensity surfaces have triaxial shape in a CDM cluster. For 
example, particles on the long axis of a prolate halo can be 
linked together, since their density is higher than the spheri- 
cal average. The subhalo catalogues we analyse in this paper 
are generated in two steps: First we use SKID with / — 4eo, 
this gives a complete catalogue of all the subhalo centers and 
also the correct subhalo properties for the smaller objects. 
Then we run AdFoF with A = 5 and combine the result- 
ing substructure catalogue with the SKID output to obtain 
the correct subhalo properties also for the larger objects: if 
AdFoF found a subhalo at the same position as SKID, the 
properties from the catalogue where this halo has a larger 
bound mass are used. The mass fraction bound to subhalos 
with N ^ 32 (the cluster centre is not considered a subhalo) 
is given in Table [3 Using the AdFoF or the SKID I = 4eo 
catalogue alone gives about 20 percent smaller values. Us- 
ing SKID with / = 1.5eo underestimates the masses of the 
biggest subhalos which dominate the bound mass fraction, 
and the results are as much as a factor of two below the 
quoted values. 

To check for systematic errors in the substructure cat- 
alogue constructed in this way, we confirmed that the sub- 
structure mass function and the number density profile of 
one cluster {D9) remains the same when we construct the 
substructure catalogue in two alternative ways: The first 
alternative catalogue was constructed by co mbining three 
SKID outputs with I — 1.5,4 and lOeo as in iGhigna et alJ 
i2000tl . the second alternative was the combination of two 



SKID outputs with / = 1.5 and 4eo and a one A = 5 Ad- 
FoF output. We found that the / = 1.5eo SKID does not 
find additional structure, the I = 4eo contains all the small 
subhalos down to the minimum number of 10. By compar- 
ing the final halo catalogue of cluster D12 to regions of the 
density map of this cluster (Figure we checked that no 
subhalos were missed and that no non-existent halos were 
included. 



3 CLUSTER SUBSTRUCTURE 

We identified subhalos within the virial radii of our six clus- 
ters at redshift zero, the algorithms used are described in 
section 12.31 At the highest resolution we found over 5000 
subhalos (^32 particles) inside the virial radius of the most 
massive cluster. 

3.1 Spatial Antibias and Convergence Tests 

In this section we study the convergence of substructure 
properties, including density profiles, cumulative mass func- 
tions and relative number density profiles (Figure First 
we check if these properties change with varying force and 
time resolution, i.e we compare D6 and D6h; D9 and D9lt; 
F9, F9cm and F9ft. The only slight difference we found 
is in the relative number density profile: the better force 
resolution in D6h leads to a few more surviving substruc- 
tures near the center (4 subhalos within 10 percent of the 
virial radius), run D6 has no subhalos within the same ra- 
dius. Therefore the original numerical overmerging problem 
iMoore. Katz fc Lakel[l99^ due to insufficient force reso- 
lution is not the limiting factor anymore, except near the 
center of the halos (r < O.lrviriai). 
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Figure 1. Density map for run D12 out to tlie virial radius. This cluster is prolate with a 3:1 inajoriminor axis 
ratio. Higher resolution color pictures and a mpeg movie of the formation of cluster C9 can be downloaded from 
: / /www-theorie.physik.unizh.ch/ ~diemand / cl usters / 1 



The amount of substructure that our simulations can 
resolve is mostly limited by mass resolution. Subhalos have 
very high phase space densities, i.e. relatively short relax- 
ation times. Numerical two body relaxation due to finite 
mass resolution he ats up their cores and makes them less 
dense jPiemand et al. 20048) . The difference in central den- 
sity is about a factor of two between subhalos resolved with 
500 and 4000 particles (see panel (a) of Figure Inland also 
iKaz antzidis ct al. (200"!) where subhalo profiles from clus- 
ters DQ, D9 and D12 and their evolution are presented). 
Subhalos with even less particles show this effect more 
strongly and have much shallower density profiles. These 
are l ess resistant against tidal stripping and total disruption 
l|Moorc ct al. 1996). 

Figure|21shows substructure properties of the same clus- 
ter, D, simulated at different mass resolutions with A'viriai 
= 205k, 1.7M, 6M and 14M. Panels (c) and (d) of Figure 



121 show the cumulative mass function including all subha- 
los with more than 10 particles. Resolution clearly affects 
the numbers of subhalos at the limiting mass of 10 par- 
ticle masses {rrip), however the amount of surviving sub- 
structure converges at a mass of about 100 nip for the D6h 
run. In analogy with t he convergence in density profiles (see 
IPiemand et aljf2004bl and references therein) we do not ex- 
pect that this number is valid for a large range of mass 
resolutions and it is possible that the high resolution mass 
functions are only complete above a mass of a few hundred 
particle masses, especially in the inner region. We usually 
include all subhalos with at least 32 bound particles for the 
analysis presented in this paper, and we will always show 
how the results depend on this minimal number of particles 
(in most cases the inffuence is small). 

Panel (b) of Figure |2| shows the number density of sub- 
halos in spherical bins relative to the number density within 
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the virial radius < n^i^iai >= A'^sut/V^iriai- The first bin is 
centered on the innermost subhalo (the cluster center is not 
considered as a subhalo), so the first data point also gives 
the radius rmin of the subhalo closest to the center. The size 
of each bin is set to rmin, so the first bin starts at rmin/2 
and ends at l.Srmin- Tidal disruption is most effective near 
the cluster center which leads to an antibias in the density 
profile of substructure relative to the smooth background. 
This implies that if galaxies are associated with the subha- 
los, they do not trace the matter distribution of a cluster. 
Is this antibias real or just an effect of finite resolution? 
Runs D&h, D9 and D12 have very similar relative number 
density profiles. If one only considers groups above the 10 
particle limit of D6h (i. e. above 80 nip in run D12), run 
D12 resolves about twice as many halos as D6h (920 against 
582, at the vertical line in Panel c) and it is interesting to 
see where these halos lie. They are not significantly more 
centrally concentrated, they have a very similar radial dis- 
tribution as the halos that survived in run D6h. Even the 
subhalo distribution of all subhalos in D3h (N^ 10) is very 
similar to the one of the subhalos in D12 in the same mass 
range (N ^ 640) which are resolved with 64 times more par- 
ticles. If the convergence scale depends only mildly on N, 
for example r converged PC N~^^^ as i n the case of the density 
profiles (see iDiemand et alJ l2004bl and references therein) , 
the wide range of resolutions presented here gives for the 
first time a robust confirmation of convergence in the radial 
distribution of subhalos. So the antibias in number density 
does not depend on the numerical resolution, but the higher 
resolution runs allow to measure the number density profiles 
closer to the center. 

The relative number density of subhalos can be approx- 
imated by an isothermal profile with a core shown by the 
thin solid line in panel (b) of Figure |5] 

n{r) = 2nH{l + {r/rHf)~\ (1) 

where nn is the relative number density at a subhalos scale 
radius rn . The average core radius of the distribution of clus- 
ter subhalos is rn ^ 0.37 rviriai — 2/3 rvc.max, where rvc.max 
is the radius where the circular velocity has its maximum, 
see Table ((TJ. 



3.2 Substructure abundance 

Figure 121 shows the cumulative substructure mass functions 
and inner mass functions of the six clusters which are all well 
approximated by a simple power law m~^. Here we include 
subhalos with a minimum of 32 particles, we found in the last 
section that the subhalo catalogues are complete only above 
a mass corresponding to about 100 particles. The apparent 
flattening of the slope towards this mass is due to finite 
resolution and does not indicate a shallower power law at 
lower masses. This can also be seen from the fact that around 
m = 10~^Mviriai the larger halos and run _D12 (i.e. those 
with better relative mass resolution) have steeper slopes. If 
hierarchical merging should produce subhalo mass functions 
th at do not dep e nd on the mass of the parent halo (as shown 
in iMoore et alJ (|^9^, see also Section 13.51 of this paper) 
the natural outcome is an power law: If one simply 

adds two equal halos the amount of substructure above any 
flxed absolute mass doubles, the remnant has now twice the 



mass and it only has the same amount of substructure at a 
flxed relative mass if the mass function of the progenitors 
was . The mass function of isolate d field halos is also 
close to a power law of slope m" -1 fe.g. lJenkins et alJl200ll : 
iReed et ahlboOSTl . Thus tidal stripping acts to lose mass in 
such a way that the overall mass function slope does not 
change. The conspiracy is such that stripped halos move 
down the M versus Hcmax plane such that they follow the 
line for fleld halos fGhiena et al. 1998). 

The cumulative substructure velocity functions (see 
Figure|3 Panel (c)) gives the number of subhalos with max- 
imum circular velocities above a given value. The virial the- 
orem Wc,max oc Mhaio/i? oc -^haio/Mj^^j^ leads to a simple 
scaling Afhaio oc fc^max for fleld halos. This relation is also 
a good approximation for subhalos, ev en if they lost most 
of their mass due t o tidal stripping (jOhigna et al.l Il998l : 
iKravtsov et alj|2004ll . Since the cumulative mass function 
goes like , we expect the cumulative mass functions to 
follow a UtTmjjx power law. This is true in a wide range of ve- 
locities. Towards the resolution limit the velocity functions 
also become shallower, but this is due to the same numeri- 
cal effect as in the case of the mass functions. The scatter 
in the substructure abundance is large at the high mass end 
(a factor three) where the mass functions depend on a small 
number of massive objects. At intermediate and small sub- 
halo masses (< 10~*M„ir) the scatter is within a factor of 
1.7. 



3.3 Subhalo velocity distribution 

3.3.1 Velocity bias 

Figure 13 panel (d) shows the 3D velocity dispersion of the 
smooth particle background and subhalos. We measured the 
dispersion profile for each individual cluster, then we aver- 
aged the values in each bin over all six clusters. The subhalos 
dispersions are not weighted by mass, each subhalo has equal 
weight. In a radial range from 0.1 rviriai to 0.4 rviriai the sub- 
structure halos have a higher 3D velocity dispersion than the 
background: b — asub/cTDM is b = 1.25 ± 0.08. The velocity 
bias of all subhalos within the virial radius b — ag^t/o'DM is 
b = 1.11 ± 0.04. The plotted and quoted errors are the scat- 
ter in our sample of six clusters and they are much larger 
than the Poisson noise in the estimated values of asub- 

A negative veloc ity bias was first considered by 

ICarlberg fc CouchmanI lll989ll as a possible way of reconcil- 
ing low cluster masses with a high matter de nsity universe. 
Hints for positi ve bias (fe > 1) were found bvlGhigna et alJ 
Jl998l) and also lCoUn. Klvpin. fc Kravtsovl tQQ(t\ who com- 
bined 12 clusters containing 33 - 246 resolved subhalos to 
obtain a sufficiently large subhalo sample. The first simu- 
lation with sufficient resolution (about 5 million particles 
within the virial radius) to construct a reliable subhalo ve- 
locity dispersion profile from one object was analysed in 
[ChiEna ct al. (2000). They found b = 1.2 - 1.3 in their in- 
nermost bin, which goes from to 0.25 rviriai, and a small 
(< 1.10) positive bias for the entire cluster. 

The bias is independent of subhalo mass, for example 
including only halos above 5 x lO'^Mviriai (979 subhalos or 
about 8 percent of the subhalos with N ^ 32) also gives 
b = 1.11 ± 0.04. And for halos above 10"*Mviriai (only 474 
halos or 4 percent) b — 1.10 ± 0.05. The velocity bias does 
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Figure 2. Substructure properties at different mass resolutions: (a) Average density profiles of 5 subhalos with masses close to 2.9 X 
10~*Mvirial resolved with about 4,000 particles in run D12 and 500 in run D6h. (b) Relative number density of subhalos with different 
mass and force resolution, fitted by an isothermal profile with a core Q The thick line is the density profile of the DM particles, (c) 
Cumulative mass functions of substructure within rvirial including halos down to 10 nip. (d) Inner cumulative mass functions, same as 
(c) but only including halos within 0.5 r^^^^^i. 



not depend on resolution: In the radial range from 0.1 rviriai 
to 0.4 rvirial the values lie within h — 1.16 and b = 1.25 for 
all simulations of cluster D and there is no clear trend with 
resolution. 



3.3.2 Anisotropy of subhalo velocities 

In the radial and tangential velocity dispersions the bias is 
very similar as in the three dimensional dispersion. This can 
also be seen from the anisotropy parameter /3 = l—0.5a^ /a^, 
(Panel (d) in Figure l^l: The anisotropy is very similar for 
subhalos and background particles, only in the inner re- 
gion the subhalo velocities are slightly more isotropic than 
those of the particle background. From r = to rvirial the 
anisotropy /3 grows roughly linear with radius: /3 ~ 0.35r. 



For the average particle anisotropy /3 ~ 0.35r^/'' seems to 
fit the data better. 



3.3.3 Subhalo dynamics 

Here we investigate if the spatial and velocity distribution 
can be a steady-state solution of the coUisionless Boltz- 
mann equation (CBE) or if a supply of infalling structures 
is needed to maintain the state of the system observed at 
z — 0. We neglect the small anisotropy and assume spheri- 
cal symmetry, then the int egral of the second momen t of the 
CBE, the Jeans Equation jBinnev &i Tremainelll98^ . reads 

Psub(r)crr,sub(7-) = / Psub(r) 5 — dr (2) 
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Figure 3. Substructure properties of the six clusters. Only halos with at least 32 bound particles are considered, (a) Cumulative mass 
functions of substructure within Tviriai- (b) Inner cumulative mass functions, including halos within 0.5 rvirial- (c) Cumulative number of 
subhalos as a function of their circular velocity, (d) 3D velocity dispersion of sub halos (circles) and dark matter background (squares) 
as a function of radius. Averages over all 6 cluster profiles, normalised to the maximum circular velocity. Error bars show the scatter 
between the clusters. Poisson errors due to small number of subhalos per bin are smaller than 0.05 and are not included. The average 
of the anisotropy parameter /3 = 1 — 0.5a^/a'^ is also plotted for the subhalos (circles) and the particles (squares). The particles are on 
slightly more radial orbits than the subhalos. The dotted lines are fitting functions, see text for details. 



where c gives the size of the system, p^uh and (Jr_sub are the 
density and the one- dimensional dispersion of the subhalos 
and M(r) is the cumulative total mass. A similar equation 
for the dark matter background is obtained by using density 
and dispersion of the dark matter instead. 

The six clusters can be approximated as NFW profiles 
jNavarro et al.lll99ef) with a mean concentration of about 
cnfw = 7 (see iDiemand et al.ll2004bl) . Using this average 
dark matter density profile the ct^dmC'') S'om Equation J^J 
fit the measured values (Figure very well. For the radial 
density profile of the subhalos we use Equation Q, with 



ru 



2/3 



the mean of 



is about 0.57 



bth = 



ctt, suhjr) 

0"r,DM('r) 



PDM(r) Psui,{r)—p^dr 



-I 1/2 



(3) 



The expected bias is 



We use a cut off at c = 2 r-viriai, at this radius the slopes 
of psub and Pdm become similar and the bias should vanish. 
Figure |1| shows the predicted and measured velocity bias 
and simple power law fit to the measured average velocity 
bias: bat ~ 1-12 x (r/rviriai)"'''^- feth is very close to the 
measured velocity bias, just in the inner region bth is too 
large. This means that the subhalo-background system is 
close to a steady-state equilibrium configuration. 

Therefore we expect the non-equilibrium processes to 
be subdominant. The net infall of subhalos can be quanti- 
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Figure 4. Velocity bias profile. Circles give the average bias of 
the six clusters. The dashed line is the bias calculated from the 
Jeans equation using the different density profiles of subhalos 
and background particles and assuming that the two are in dy- 
namical equilibrium. The dashed-dotted line gives a simple power 
law fit DC r~^'^ to the average bias. 



fied from the asymmetry of the radial velocity distribution 
of subhalos near the virial radius: The distributions are sym- 
metric in the inner and outer part of the clusters and there 
is no net infall of subhalos at z=0. Another non-equilibrium 
process is the disruption of subhalos. The fraction of subha- 
los that are disrupted is small (see also Section . about 
0.02 Gyr"^ for subhalos with N > 100. In the inner 40 per- 
cent of the halo the fraction is bigger, about 0.13 Gyr~^. 
This could be the reason why the steady-state solution over- 
predicts the velocity bias near the center. 



3.3.4 Higher moments of the velocity distribution 

In the last subsection we found that the second moment of 
the subhalo velocity distribution is consistent with a steady- 
state solution, where the subhalos have a spatial antibias. 
Now we consider the next higher moments of the velocity 
distributions of subhalos and particle background. In the 
radial range where the velocity bias is large (0.1 rviriai to 0.4 
^virial) the shapcs of these velocity distributions are very dif- 
ferent (Figure |SJ. There are many less subhalos with small 
velocities (top panel), also the fraction of subhalos with 
low velocity components is smaller for the particles (bot- 
tom panel). While the particle velocity distribution is close 
to a Maxwellian, this is not true for the subhalos. The sub- 
halo velocity histogram is flat-topped, it has smaller fourth 
moment than the Maxwell distribution, i.e. a negative kur- 
tosis k =< V* > / < >^ — 3 = —0.7. We also calculated 
the first two non-t rivial, even^ Gauss-Hermite moments /i4, 
h6 jGerharJ l993?^. In this radial range (0.1 - 0.4 rviriai) we 
get /i4 = —0.068 and h^ — 0.0013. The advantage of Gauss- 
Hermite moments over simple higher order moments is that 




1 1.5 

V / max V . 




The odd moments are zero for symmetric functions. 



V / max V . 

circ 

Figure 5. Velocity distribution of inner subhalos (circles) and 
particles (squares). Average of six distributions from different 
clusters. Subhalos and particles between r = O-lr^iriai and 
0.4rviriai are included. Velocities are normalized to the maximum 
circular velocity Vc,max of each cluster. Solid and dashed lines are 
Maxwellian distributions with the correct second moment. 



they are not very sensitive to the wings of the distribution. In 
galaxy clusters these outer parts of the distribution are hard 
to de termine exactly due to interlopers ivan der Marel et alJ 
l200nh . 

In Figure |S| we plot the velocity histogram further out 
(0.5 rviriai to rviriai). Now the second moments of the particle 
and subhalo velocities are much closer {b = 1.10), but the 
shapes of the velocity distributions of subhalos and particles 
are stiU different: k = -0.60, = -0.031 and he = -0.025. 
For all subhalos within rviriai we find b — 1.11, k = —0.48, 
h4 = -0.034 and he = -0.012. 

Both the inner (Figure|3l and outer (Figure|S|l subhalos 
show an excess of high-velocity substructures between Hcmax 
and 1.5 Uc.max . Many of these high-velocity subhalos are on 
very radial orbits. When we exclude subhalos with absolute 
values of the radial velocity component larger than Wc.max 
the excess disappears and the speed distribution follows the 
Maxwellian distribution of the background particles above 
Wc,max . The large fraction of subhalos with very high radial 
velocities is also evident in the radial velocity distribution 
(not shown): both in the inner and outer part of the clusters 
the distribution has a very negative kurtosis of fc = —0.9. 
Also note that the radial velocity distributions are symmet- 
ric, there is not net infall of subhalos at z=0. 

The shape parameters depend weakly on the lower mass 
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Figure 6. Same as Figure|S] but including subhalos and particles 
between r = O.Srvirial and Tvirial- 



threshold, including subhalos above 5 x 10 ^Afviriai instead 
of 32 mp yields: b = 1.11, k = -0.44, M = -0.016 and 
ha — —0.022. There are 979 subhalos above this threshold 
in our six clusters, which is only 979/12027 = 0.039 of the 
N ^ 32 subhalo sample, but this is still enough to determine 
the shape of the velocity distribution. All of these subhalos 



have bound masses of more than 1.2 x 10'"M, 



©• 



3.4 The origin of tlie subhalo biases 

The physical mechanism that generates the differences in 
the spatial and velocity distributions of particles and sub- 
halos is most likely the tidal destruction of subhalos in dense 
environments. The efhcency of tidal stripping and tidal dis- 



on the orbital enerf 


lY of the subha- 


iTaffoni et al.ll2003l: 


iKravtsov et alJ 



120041) . Therefore it offers a natural explanation for the lack 
of slow subhalos; at a fixed radius the orbital energy is pro- 
portional to the square of the velocity and tidal disruption 
could remove a large fraction of the slow subhalos producing 
a distribution like the one given in the top panel of Figure 
IS] 

The tidal disruption of subhalos must occur very early 
in the evolution of the cluster. iGhigna et aL (i20OTi) are able 
to identify the remnants of 60 to 70 percent of all cluster 
progenitor halos (A'^ 100 at z=3) with subhalos at z=0. 



. From the halos identified at z = 1 an even larger fraction 
survives (more than 80 percent). For run D6h we performed 
the same test and get very similar numbers. We link pro- 
genitor halos with a halo at z = if at least four particles of 
the progenitor are bound to the subhalo at z — and find 
descendents for 83 percent of the progenitor halos identified 
at z—2. 

However a significant fraction of subhalos may have 
been destroyed prior to this epoch. From the halos with 
N ^ 100 identified in the high resolution region of run D6h 
at z=7.2 and z=4.3, we can associate only about 60 percent 
with z — subhalos. At this early stage tidal disruption 
seems to act as a physical selection process which allows 
only halos with high enough orbital energies to survive as 
todays subhalos. This causes the spatial antibias and the 
positive velocity bias of substructure. 

Note that it is important to have a larger minimum 
number of bound particles in the early subhalo sample {N ^ 
100) than in the final subhalo catalogue (N ^ 10) if one 
wants to quantify disruption: If we would use the same N 
at both times then we would get a much higher 'disruption 
rate', but we would mostly measure the amount of subhalos 
that were tidally stripped below this threshold number of 
bound particles but not necessarily disrupted. This caveat 
would have a big infiuence since about half of the considered 
subhalos have a bound mass between Nm^ and 2Nmr,. 



3.5 Comparison with galaxy size halos 

The four galaxies in our sample are resolved with 1.7 to 3.8 
million particles, so the relative mass resolution is lower than 
for the clusters. However, there is enough resolved substruc- 
ture to compare its abundance and the radial distribution to 
the results from the cluster runs. We make the comparison 
with cluster D6h which has similar relative mass and force 
resolution as the galaxies. We also give the results for the 
same cluster with eight times better mass resolution (run 
D12) to get an impression how the results might change if 
we also had higher resolution for the galaxy halos. There 
may be a hint that the galaxies have slightly less substruc- 
ture than the clusters, but we need to increase the resolution 
in the galaxy simulations in order to verify this result. 

Galaxy G2 had a recent major merger at z ~ 0.2, at z = 
this merger is finished, the core has no more visible signs of 
dynamical activity. The concentration of this galaxy is lower 
Cvc,max = rviriai/rvcmax ^ 3.6, probably due to the later 
formation in this recent merger. The other three galaxies 
had no more major mergers since at least z ~ 0.2 and their 
Cvc max are between 5 and 6.5. 



3.5.1 Substructure abundance 

Despite the fact that clusters form much later than galax- 
ies in hierarchical structur e formation, they ha ve very sim- 
ilar subhalo mass function. iMoore et al.l 1^9^ showed this 



^ The fraction of subhalos that merge with the central object (i.e., 
end up within an assumed radius of about 0.015 rviriai) are always 
below 5 percent and can be neglected in this context. But it is an 
impo rtant fraction if one considers the most massive progenitors 
only jGhigna et ai]l2000h . 
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by comparing two SCDM halos. iDe Lucia et alJ i2004) con- 
firmed this recently for several ACDM halos, but at a reso- 
lution of less than a million particles inside the virial radius. 
Figure |7| shows the subhalo abundance in the four galaxies 
and in the cluster D. The velocity functions (Panel a), mass 
functions (Panel c) and inner mass functions (Panel d) are 
all quite close to those of the reference cluster run DQh. 

The substructure abundance is largest in galaxy G2, it is 
as high as in run D&h. This halo formed recently in a major 
merger at z ~ 0.2, which is a typical formation history for 
cluster size halos rather than for galaxies. The other three 
galaxies have about 30 percent less substructure than G2 
and D&h. Therefore the amount of substructure depends 
weakly on the mass of the parent halo, but the difference 
appears to be comparable to the scatter within parent halos 
of a fixed mass. 



3.5.2 Radial distribution 

The relative number density profiles (Panel (b) of Figure 
of the galaxy subhalos ar e more centra llv concentrated than 
those of cluster subhalos Lucia et al...2004l. Sm aller ha- 
los have higher concentrations |Navam^^Llll99^ and are 
therefore more resistant against tidal disruption. However 
the subhalo number density also shows a clear antibias with 
respect to the dark matter density. 

The density profile that fits the cluster subhalos distri- 
bution (Equation is a good approximation also for the 
galaxy subhalo number density profile. Now the core radius 
is a smaller fraction of the virial radius (rn — 0.14 rviriai) 
because galaxy subhalos are more centrally concentrated. 
Note that rn is again about two thirds of the radius where 
the circular velocity is maximal, this is the same fraction as 
for the cluster subhalos. Therefore scaled to rvc.max galaxy 
and cluster subhalos number density profiles are the same. 
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Figure 8. Projected relative number surface density profile of 
subhalos averaged over the six clusters: Circles include all (12'023) 



and plus signs only halos with i;c,max > O.OQfc, max, main - The core 
of the main halo, the 'cD galaxy', is always included in the first 
bin. The projected dark matter density is plotte d with squares. 
Crosses are the data from the CNOC surve y {Ca rlbcrg ct al, 
,^2^, stars are the Coma cluster data from iLokas MamoE 
i2003r . We normalise the curves so they match at rviriai. 



this velocity. But the number of Milky Way sate llites with 
«c,max > 0.05t;c,max,parcnt Is Only 10 ijMoore et alllQQft 'l. Var- 
ious soluti ons to this issue have been proposed in the liter- 
ature (e.g. lStoehr et. aJl20oi iKravtsov et al.ll2004l . 



4 COMPARISON WITH OBSERVATIONS 
4.1 Substructure abundance 

iDesai et al.l l|200^ measured galaxy circular velocity func- 
tion in 34 low-redshift clusters and found that these func- 
tions can be approximated by a power-law oc v'^^'^^^. In 
CDM cluster simulations they found a logarithmic slope of 
— 3.4±0.8. Our higher resolution simulations show that these 
slopes are rather on the steep side of the given range. Figure 
^ shows that the cumulative velocity function has a slope 
of about —3, where we expect the sample to be complete. 
For the differential circular velocity function this gives a 
slope of —4, which is not consistent with the observed slope 
of —2.5. Accounting for the effects of the baryons could 
reconcile CDM simu latio ns with the observ ations, see e.g. 
ISprineel et"ai] (|20o3) and lDesai et alj ll2004li . Realistic gas- 
dynamical cluster simulations will eventually resolve this is- 
sue. 

The same problem is more severe when the host halo 
is a galaxy and not a cluster. The steep circular velocity 
function of CDM halos predicts over 100 subhalos with cir- 
cular velocities above 5 percent of the parent halo circular 
velocity, i.e. above 10 km/s for a Milky Way size halo. Our 
highest resolution cluster ^49 has over 300 subhalos above 



4.2 Spatial distribution 

For comparison with observed spatial and velocity distribu- 
tions of galaxies in clusters we 'observe' the six simulated 
clusters along three different line of sights (LOS) (the x, y 
and z axis) and average over these LOS. We then take the 
sample averages to get mean values and an estimate of the 
scatter. The results are shown in Figure Eland El 

The number surface density is plotted at the midpoints 
of equal bins in projected distance from the densest region 
of the cluster. The innermost bin starts at J? = and there- 
fore always one additional subhalo, i.e. the core is counted 
as the cD galaxy of the cluster. The projected number den- 
sity is flat near the center, just like the 3D number den- 
sity in Figure 121 The total sample contains 12'027 subhalos 
with at least 32 bound particles from the 6 high resolution 
clusters. In the Coma cluster a number density profile for 
a comparable number of galaxies (985) can be measured 
(Lokas & Mamon 2003), this profile (plotted with stars in 
Figure [SJ is steeper than the subhalo profile and follows 
rather closely the expected d ark matter profile of a CDM 
cluster. ICarlberg et al.l lll997li give the surface density pro- 
file of a sample of galaxies combined from 14 clusters ob- 
served in the CNOC cluster survey. The sample contains 
1150 galaxies, including background and goes out to 2r2oo, 
i.e., per cluster there are about 50 galaxies. Therefore this 
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Figure 7. Substructure properties of four galaxy halos: Panel a): Cumulative number of subhalos as a function of their circular velocity. 
Panel b): Relative number density of subhalos and of all DM particles (see Section l3.1l for details). Panel c): Cumulative mass functions 
of substructure within Tviriai. Panel d): Inner cumulative mass functions, including halos within 0.5 rvirial- AH halos with at least 10 
bound particles are included in these plots. The solid lines show the four galaxies, the dashed line is a cluster halo at similar resolution 
and the dashed dotted line is the same cluster at eight times higher force resolution for comparison. 



magnitude limited sample should be comparable to the most 
massive 300 subhalos in our sample. 

We selected the subhalos with m > 2 x Alyiriai and 
get a sample of 238 halos, their surface density profile is plot- 
ted with triangles in Figure |5| The profile does not change 
much, just in the innermost bin the values rise, due to the 
relative importance of the 'cD galaxy'. Selecting subhalos 
by peak circular velocity Ucmax > 0.09 i;c, max, parent gives a 
sample of 291 halos with a similar surface density profile. 

The observed n umb er surface density pro files from 
lll997[l and lLokas fc MamonI (l2003l) (and also 



C^lbe^^^l 



Beers fc Tonrvl lll986ll and lMerrifield fc KenI ^9^)) are sig 



nificantly steeper than in the CDM clusters. To correct the 
subhalo number density in the inner four bins upwards to 
match the observed values one needs to add a number of 



subhalos similar to the total number within the virial radius 
of each cluster, but preferentially more subhalos closer to 
the cluster center. We discuss the implications of this result 
in the conclusions. 



4.3 Subhalo velocities 

The velocity bias b ~ 1.12 ± 0.04 would lead to dynamical 
cluster mass estimates that are about 20 percent too high if 
cluster galaxies reside in CDM subhalos. By comparing with 
cluster mass estimates from gravitational lensing it could be 
noted the dynamical estimated are too high, but it is very 
difficult to obtain estimates with small enough uncertainties 
with both methods. S uch a comparison was performed by 
ICvpriano et al.l i200 J) , finding that dynamical masses are 
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indeed biased by 1.20 ± 0.13 in a sample of 14 clusters, but 
the effect only comes from the massive clusters ((t„ > 1122 
km/s), which show large mass differences 1.54 ± 0.19, while 
the smaller clusters show no bias. 

FigureOshows the projected moments of the CDM sub- 
halo velocity distributions and the inner and outer distribu- 
tion of line of sight velocities averaged. The ploted values 
are averages over the six cluster halos and over three differ- 
ent projections. The velocity moments for the dark matter 
background are also plotted for comparison, a s i milar anal- 
ysis was presented bv lSanchis. Lokas fc MamonI (|20o3)- 

In contrast to the spatial distribution the velocity dis- 
tribution of CDM subhalos agrees surprisingly well with 
current observations of cluster galaxies. In the grand to- 
tal velocity distribution of the CNQC survey a negative 
hi = -0.015 ± 0.005 was found Jvan der Marel et alj|200oll 
and /i6 = -0.028 ± 0.006. We get k = -0.44, M = -0.016 
and /i6 = —0.022 using all subhalos with bound mass larger 
than 5 x 10~^A/viriai. There are 1152 subhalos above this 
threshold in our six clusters. This agreement between simu- 
lations and observations may be fortuitous since the spatial 
distribution of galaxies is different and probably due to de- 
struction of low energy central subhalos. Also, in the Coma 
cluster the velocity distribution seems to be more flat topped 
compared to a Gaussian: The kurtosis is negative in most 
radi al bins, the value s scatt er around k ~ —0.5 (see Figure 
3 in lLokas fc MamonI (|2D03)). The uncertainties in the mea- 
surement of velocity moment profiles are still quite large and 
a comparison with the projected moments from Figure |S| of 
the CDM subhalo velocities is not feasible yet. 



5 CONCLUSIONS 

We analyse the substructure within six very high resolution 
cold dark matter simulations of galaxy clusters and four sim- 
ulations of galaxies. We have addressed several open issues 
raised in the introduction regarding the results of high reso- 
lution simulations of individual halos within the concordance 
CDM model. Our conclusions can be summarised as follows: 

(i) The spatial distribution of subhalos in cold dark mat- 
ter simulations of galaxies and clusters is antibiased with 
respect to the mass. Although this behavior was found by 
other groups, we demonstrate that this result is robust and 
does not change as we increase the resolution. We show that 
this antibias most likely results from a population of early 
halos that are tidally destroyed in the dense protocluster en- 
vironment and within the central regions of the final cluster. 

(ii) The surviving population of subhalos have a posi- 
tive velocity bias that increases towards the center of the 
halos. The subhalo velocity distribution functions are non- 
Gaussian, they are 'flat topped', especially in the inner re- 
gion: There the kurtosis is fc = —0.7 and the fourth Gauss- 
Hermite coefficient hi — —0.068. 

(iii) The spatial anti-bias and the positive velocity bias 
of the subhalos are consistent with a steady-state solution 
of the Jeans equation. Subhalos are a hot, more extended 
component in equilibrium with the potential generated by 
the smooth particle background. 

(iv) The mass and circular velocity distributions of sub- 
halos in our highest resolution simulation show the same 
power law slopes as in lower resolution versions, but are 



steeper at the low mass end. It is not clear that convergence 
in the number of subhalos has been reached below a scale of 
a few hundred particle masses. 

(v) Cluster and galaxy mass halos simulated at the same 
resolution have similar substructure abundances. The scat- 
ter in the circular velocity and mass functions is a factor 
of three at the high mass end, but falls to just 1.7 at lower 
masses. 

(vi) An observational comparison with CNOC cluster 
data and the Coma cluster shows that the galaxy popu- 
lation traces the smooth dark matter background, but not 
the predicted halo population. This is most likely due to 
overmerging in the central region of the simulations and we 
are probably missing a factor of two in the subhalo popu- 
lation. The baryonic cores of these disrupted subhalos may 
survive intact if dissipational processes increase their densi- 
ties sufficiently. Also a greatly truncated dark matter halo 
may survive in this case. 

This latter statement is the most profound conclusion 
of this work. The spatial distribution of cluster galaxies is 
significantly different from the distribution of subhalos in 
dark matter simulations. Either the model is incorrect or we 
have reached a fundamental limit to this type of pure dark 
matter simulation. Here we explore the latter possibility and 
the implications for the morphology density relation. 

It is likely that disk galaxies do not significantly mod- 
ify the overall potential provided by the baryons and dark 
matter. Whereas a disk-disk merger would funnel gas to the 
central region, forming an elliptical galaxy with a signifi- 
cantly deeper potential and a effective rotation curve that 
is at least isothermal, or possibly Keplerian in the center 
iRomanowskv et al.ll20o3) . Thus we expect that an ellipti- 
cal galaxy would most likely survive at any position within 
the cluster, albeit with a greatly truncated dark matter halo. 
Late type spiral galaxies are unlikely to survive within the 
central regions of clusters (or their progenitors) and will be- 
come physically overmerged to form the cD halo of diffuse 
light. 

If the CDM paradigm is correct then we are missing 
close to a factor of two of the 'galaxy' population as asso- 
ciated with subhalos, increasing to a factor of five within 
the inner 10% of the cluster. It is possible that simulations 
with more than 10^ particles per system may resolve more 
central subhalos and calculations this large will be possi- 
ble in the future. In this case, the velocity bias should de- 
crease as we resolve more halos/galaxies in the central re- 
gions. However, from our convergence study we find very 
few new halos in the central cluster regions as we increase 
the resolution by a factor of ten. This implies that we have 
reached a physical limit to DM-only simulations and that 
any loss of subh alos in current simulations is due to ph ysi- 
cal overmerging JWhite fc ReeJl978l:lMoore et al]ll996l) . In 
this case progress in this area can only be made by including 
a realistic treatment of hydrodynamics and star-formation 
such that realistic disks and elliptical galaxies can be fol- 
lowed within the appropriate cosmological context. 

The survival or disruption of a galaxy depends on an 
intricate balance between the progenitors dark halo struc- 
ture and the effects of dissipation. Sa-Sb galaxies must lie on 
the borderline between survival and disruption in the cluster 
environment. The morphology-density relation may simply 
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Figure 9. Panel a): Average line of sight velocity dispersion of subhalos and particles as a function of projected distance from the center. 
Panel b) and panel d): Average line of sight velocity distributions of subhalos and particles, for projected radii smaller (b) and larger (d) 
than 0.4 rvirial- Solid and dashed lines are Gaussians with a second moments fitted to the subhalos (solid) and to the particles (dashed). 
Fourth order Gauss-Hermite approximations to the subhalo velocity distribution functions are given with dashed-dotted lines. Panel c): 
Average kurtosis (with error bars) and fourth Gauss-Hermite moment (without error bars and multiplied by a factor of ten for clarity) 
of the line of sight velocity components of subhalos and particles as a function of projected distance from the center. The error bars in 
panels a) and c) give the scatter within the six clusters. 



reflect the fact that the disks are preferentially destroyed in 
the central regions of clusters. However if the CDM model 
is correct one needs to preferentially form ellipticals in high 
density regions before the cluster forms. The fact that the 
observed galaxy distribution follows the dark matter distri- 
bution implies that no overmerging of galaxies has taken 
place. It is insufficient to take disks and destroy them in 
the cluster cores since this would give rise to a cored galaxy 
distribution. 

The fact that 40% of halos identified at z—7 can not 
be associated with a subhalo at z—Q, or have not merged 
with the central cD, implies that they have merged into the 
smooth particle background. If these objects can be associ- 
ated with surviving galaxies, it implies a strong age-radius 



dependence for galaxies within clusters. At the cluster cen- 
tres over 80% of the galaxies must have formed prior to z=7. 
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